COBE Preprint # 96-06, Submitted to ApJ Letters 



2-Point Correlations in the COBE 1 DMR 
4- Year Anisotropy Maps 

G. Hinshaw 2 - 3 , A.J. Banday 2 ' 4 , C.L. Bennett 5 , K.M. Gorski 2 - 6 , A. Kogut 2 , 
C.H. Lineweaver 7 , G.F. Smoot 8 ' 9 & E.L. Wright 10 

ABSTRACT 

The 2-point temperature correlation function is evaluated from the 4-year COBE 
DMR microwave anisotropy maps. We examine the 2-point function, which is the 
Legendre transform of the angular power spectrum, and show that the data are 
statistically consistent from channel to channel and frequency to frequency. The most 
likely quadrupole normalization is computed for a scale-invariant power-law spectrum 
of CMB anisotropy, using a variety of data combinations. For a given data set, the 
normalization inferred from the 2-point data is consistent with that inferred by other 
methods. The smallest and largest normalization deduced from any data combination 
are 16.4 and 19.6 fiK respectively, with a value ~18 fiK generally preferred. 
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1. Introduction 

The detection of large angular scale anisotropics in the Cosmic Microwave Background (CMB) 
radiation was first reported by the COBE-DMR experiment in 1992 (Smoot et al. 1992; Bennett 
et al. 1992; Wright et al. 1992; Kogut et al. 1992). The initial detection was based only on the 
first year of flight data. Since that time the DMR Team processed and analyzed the first two years 
of data and found results to be consistent with the first year results (Bennett et al. 1994, Gorski 
et al. 1994, Wright et al. 1994). We have now processed and analyzed the full 4-years of DMR 
observations: this paper is one of a series describing the results of our analysis. The maps and an 
overview of the scientific results are given in Bennett et al. (1996). 

In this paper we analyze the anisotropy in the 4-year DMR maps using the 2-point correlation 
function as a measure of the angular power spectrum. The COBE-DMR experiment was designed 
to measure the CMB anisotropy on angular scales of > 7°, corresponding to spherical harmonic 
multipole moments of order £ < 30. The DMR has produced full-sky maps of the CMB temperature 
at each of three frequencies, 31.5, 53, and 90 GHz, with two independent channels, A and B, at 
each frequency. In principal, one can obtain an estimate of the CMB power spectrum from an 
anisotropy map simply by decomposing the map, T(9,(f>), into spherical harmonic components 
and averaging them to find the mean power per mode t. T(#, 4>) = J2e m a tmXimiPi 4>) with power 
spectrum af = 2 ^ +l J2m=-e \ a £m\ 2 - I n practice, however, there are a number of complications 
that arise. First, the need to apply a galactic cut to the data renders the spherical harmonics 
non-orthogonal, thereby coupling the at m coefficients and increasing their uncertainty. Moreover, 
since af is a quadratic form, any uncertainty in the a^ m (whether due to coupling, instrument 
noise, systematic effects and/or foreground sources) produces a positive bias in the estimate of 
the power spectrum. However, see Gorski et al. (1996) and Wright et al. (1996) for spherical 
harmonic-based analyses that account for these difficulties. 

An alternative to estimating the power spectrum is to evaluate its Legendre transform, the 
2-point correlation function. For a given power spectrum with multipole amplitudes Ci = {\a£ m \ 2 ), 
the predicted covariance between pairs of map pixels i and j with angular separation a-ij is 

C( aij ) = (TiTj) = 1-^2(21 + 1) Wf C e P*(cos atj ) (1) 

where is the CMB temperature in pixel i, the angled brackets denote an average over an 
ensemble of universal observers, We is the experimental window function that includes the effects 
of beam smoothing and finite pixel size, and P;(cos «y) is the Legendre polynomial of order £. We 
estimate the 2-point correlation function in our sky by evaluating the average product of all map 
temperatures with a fixed angular separation 



C(a) = ^WiWj TiTj j ^WiWj. 



(2) 
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where the sum is restricted to pixel pairs separated by an angle a, Ti is the observed 

temperature in pixel i after monopole and dipole (and optionally quadrupole) subtraction, and Wi 
is the statistical weight of pixel i. This statistic is straightforward to compute, and offers a quick 
test of the consistency of the power spectrum from map to map. 

In the approximation that the 2-point function can be treated as a multivariate Gaussian 
distribution, one can form a likelihood function with which to estimate the power spectrum 
normalization. In §2 we present the 2-point correlation data in the 4- year DMR maps and examine 
its consistency from map to map. In §3 we evaluate the Gaussian likelihood as a function of 
the mean expected quadrupole moment, Qrms-PSi under the assumption of a scale-invariant, 
power-law spectrum of anisotropy. In §4 we compare these results with those obtained by other 
methods and summarize our findings. 

2. 2-Point Correlation Data 

The 2-point correlation function, as given in Equation [2], is the average product of all pixel 
temperatures with a fixed angular separation. The data are binned into angular separation bins 
of width 2?6 with the first bin reserved for all pixel pairs such that i = j, the second bin for 
pairs with separation between 0° and 2?6, and so forth. For the present analysis we employ the 
maps pixelized in galactic coordinates and use the custom Galaxy cut described by Bennett et al. 
(1996), for which there are 3881 surviving pixels. To minimize cosmic variance we assign equal 
weight to each surviving pixel. 

The basic 2-point functions obtained from the single frequency maps are shown in Figure |]. 
We plot both the auto-correlation of the weighted sum of channels A and B at each frequency (the 
coefficients used to form the weighted average maps analyzed in this paper are given in Table 0) 
and the cross-correlation between channels A and B. For reference, we also plot, as a solid line, 
the auto-correlation of the weighted average of all six DMR channel maps. The error bar attached 
to each point represents the rms due to instrument noise, based on 2000 Monte Carlo simulations 
that include only instrument noise. The plot clearly demonstrates excellent consistency of the 
2-point correlations at 53 and 90 GHz, even in the absence of any galactic signal corrections. The 
31 GHz data exhibit a small discrepancy from the mean data that is primarily quadrupolar and is 
presumably due to residual galactic emission. 

The data are quantitatively tested for self-consistency by forming differences of the 
2-point functions and comparing them to simulations. The statistic for the test is defined as 
X 2 = (AC - (AC)) T • M _1 • (AC - (AC)) where AC is the observed difference between 2-point 
functions, with entries AC a = C^(a a ) — C <y2 \a a ) (a denotes an angular separation bin, (1) and 
(2) denote specific data selections), (AC) is the mean difference, computed from simulations 
described below, and M is the covariance matrix computed from simulations. For each realization 
in the Monte Carlo, we generate a single realization of a scale-invariant power-law sky with unit 
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normalization, and six noise maps, one per channel, with appropriate noise level and coverage 
(Bennett et al. 1996). We assume the noise is uncorrelated from pixel to pixel, based on the 
analysis of Lineweaver et al. (1994). It is then possible to generate an ensemble of simulated 
2-point functions for any desired auto- or cross-correlation function constructable from the DMR 
data. We generate such an ensemble for each of the six panels depicted in Figure |l| note that 
a given realization of the six combinations shares a common CMB signal. We compute x 2 as 
defined above for each possible difference and compare its value to the ensemble derived from the 
simulations. Note that our computation of the covariance matrix from simulations automatically 
includes bin-bin correlations in the definition of x 2 ■ I n no case does the observed value of x 2 
exceed the 5% confidence upper limit derived from the simulations, which corroborates the visual 
consistency of the data. 

The 2-point functions obtained from selected multi-frequency combinations of the data are 
shown in Figure |2[ We plot the auto-correlation of the weighted average map, the cross correlation 
of the 53 and 90 GHz maps, and the auto-correlation of two maps which have had residual, 
high-latitude galactic emission modeled and removed (Table |l] and Kogut et al. 1996a). Note 
the excellent consistency between the auto-correlation of the weighted average map, which is 
sensitive to all structure in that map, and the cross-correlation of the 53 and 90 GHz data, which 
is sensitive only to common structure in the maps. Note also that the two methods used to model 
and remove high-latitude galactic emission introduce only small changes in the 2-point data, and 
hence in the angular power spectrum. This observation, coupled with the fact that the Correlation 
and Combination model maps render very similar 2-point functions, supports the claim of Kogut 
et al. (1996a) that the free-free emission at high latitudes is 1) weak, and 2) approximately traced 
by the DIRBE 140 /jm map at 7° resolution. 



3. Quadrupole Normalization 

Given a power law model of initial Gaussian density fluctuations, P{k) oc k n , where P(k) is 
the power spectrum of density fluctuations as a function of comoving wavenumber k, it is possible 
to derive the corresponding angular power spectrum of CMB fluctuations, Ce = (|a^ m | 2 ) (Bond & 
Efstathiou 1987). The result is 

r = r r(l+(n-l)/2)r((9-n)/2) 
1 2 r(^ + (5-n)/2)r((3 + n)/2) K) 
For the scale-invariant case, n = 1, this reduces to Cg = 6C2/(i(i + 1)), which has one free 
parameter, the mean quadrupole moment C2. We customarily express the normalization in terms 
of Qrms-PS = V (5/47r)C2, the mean rms temperature fluctuation expected in the quadrupole 
component of the anisotropy. We determine the most likely quadrupole normalization, Q rms -ps, 
from the 2-point function by evaluating the Gaussian approximation to the likelihood function 

e -|AC T -M- 1 -AC 

CiQrms - PS) " Vdet(M) • (4) 
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Here AC T and AC are m-dimensional row and column vectors with entries 
AC a = C{a a ) — (C(a a )), and M = {(AC)(AC) T ) is the covariance matrix of the correlation 
function. Note that AC is defined differently than in §2 since we are comparing a single 2-point 
function to an ensemble here. The angled brackets denote averages over both measurement 
errors and over the ensemble of anisotropy fields implied by cosmic variance for a given Cg. We 
estimate the mean correlation and covariance matrix as a function of Q r ms-PS using Monte 
Carlo simulations described above. The simulations account for all important aspects of our data 
processing including monopole and dipole (and quadrupole) removal on the cut sky. Because of 
this subtraction, the bins of a given correlation function are not all independent so the covariance 
matrices derived from the simulations are formally singular. We invert these matrices to form 
X 2 using singular value decomposition, which permits an unambiguous identification of the zero 
modes that arise due to multipole subtraction. We then evaluate the logarithm of the likelihood, 
ln£ = — ^[x 2 + ln(det(M)) + const], in steps of 1 //K in Qrms-PS, spline the result to a resolution 
of 0.01 /uK, and identify the maximum. 

We test the likelihood method for accuracy by feeding the simulated 2-point functions into 
the likelihood function and solving for an ensemble of Q r ms-PS maxima. We define the bias in 
our method to be AQ = (Qmax) — Qin where {Qmax) is the mean of the recovered maxima and 
Qi n is the simulation input normalization. The resulting bias depends on the noise level in the 
data, but ranges from —0.2 to —0.4 fiK for the all cases except the 31 GHz data where it is ~ — 1 
IjlK. We correct for this bias in all reported results. The uncertainty we assign to Qrms-PS 1S t ne 
rms scatter of the ensemble Qmax which typically exceeds the rms of the Gaussian likelihood by 
about 10%. 

The corrected power spectrum normalization deduced from a variety of data combinations 
is given in Table |2[ The smallest and largest normalization deduced from any data combination 
are 16.4 and 19.6 yuK respectively, with values ~18 uK generally preferred. The normalization 
inferred from the 2-point function is now in better agreement with other determinations than was 
the case with the 2-year data. The change is due to data selection: with the 2-year data, we only 
analyzed the 53 x 90 GHz cross-correlation function; with the 4-year data we have analyzed many 
more data combinations, including the auto-correlation of a weighted average multi- frequency 
map which yields a normalization ~1.5 /iK higher than the cross correlation. The multi-frequency 
auto-correlation is more comparable to the data analyzed by other methods, and the 2-point 
analysis yields consistent results in that case. For a comparison, see Table 2 of Bennett et al. 
(1996). In general, the normalization inferred from the 4- year data is slightly less than we found 
after 2 years, in part because of the extension of the Galaxy cut. For comparison, the 31+53+90 
GHz auto-correlation with a straight 20° cut yields a best-fit normalization of Qrms-PS = 19 /xK. 
As shown in Table ||, the effects of further modeling and subtraction of galactic emission are less 
than 1 fiK in the normalization. 

While a likelihood analysis is capable of inferring the best-fit parameters for a given model, it 
does not say anything per se about the goodness of fit. For reference we have included in Table 
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H the values of x 2 at the maximum likelihood value of Q r ms-PS (with n = 1). Since the 2-point 
function is only approximately multivariate Gaussian distributed, our tabulated statistic is only 
approximately x 2 distributed. However, we have used our Monte Carlo simulations to compute 
the expected distribution of this statistic and find it to be approximately x 2 with a mean of ~ 70 
and a standard deviation of ~ 12. The values computed with the DMR data are very consistent 
with this distribution, implying that the data are well fit by a scale-invariant power spectrum. 



4. Conclusions 



Analyses of Q r ms-PS\n=i have also been reported by Banday et al. (1996), Gorski et al. 
(1996), Hinshaw et al. (1996), Kogut et al. (1996b), and Wright et al. (1996). All results lie 
between 15.5 and 19.6 /iK with most between 17.5 and 18.5 fiK. In general, all methods for 
analyzing a given data combination give consistent results, while there is modest dependence 
on data selection. Fortunately, the dependence on data selection does not exceed the statistical 
uncertainty due to cosmic variance and instrument noise. A reasonable conclusion to draw for the 
normalization of a scale invariant spectrum is Q rm s-PS\n=i = 18 ± 1.6 /xK, with the error quoted 
conservatively. 

We gratefully acknowledge the many people who made this paper possible: the NASA Office 
of Space Sciences, the COBE flight operations team, and all of those who helped process and 
analyze the data. 



- 7- 



Table 1. DMR Map Combination Coefficients' 1 



Map 


31A 


31B 


53A 


53B 


90A 


90B 


DIRBE 


Haslam 


31ws 


0.611 


0.389 




















53ws 








0.579 


0.421 














90ws 














0.382 


0.618 








53+90 








0.412 


0.299 


0.110 


0.179 








31+53+90 


0.049 


0.032 


0.378 


0.275 


0.102 


0.164 








Correlation 13 


0.049 


0.032 


0.378 


0.275 


0.102 


0.164 


3.364 


0.314 


Combination 


-0.185 


-0.117 


0.367 


0.266 


0.256 


0.413 


2.055 


-0.170 



a The maps are formed using the above coefficients according to the prescription: 
S« CdmrTdmr ~ ChTh — CdTd, where i = 31A, • • • , 90B is a channel index, C l DMR are the DMR 
map coefficients given above, Tomr are the DMR maps in fj,K of thermodynamic temperature, Ch 
is the coupling coefficient to the Haslam map, in fiK / K, given above, T# is the Haslam map, in 
K, Co is the coupling coefficient to the DIRBE 140 fim map, in ^iK / (MJy/sr), given above, and 
To is the DIRBE 140 fim map, in MJy/sr. The resulting map has units of fiK, thermodynamic. 

b The coefficients give the most sensitive combination of the 31, 53 and 90 GHz data, in 

thermodynamic units. The Haslam map is used to model synchrotron emission, it was fit to 

the DMR data under the assumption that its spectral index was (3 S = —3.0. The DIRBE 140 //m 

map is used to model both free- free and dust emission, see Kogut et al. (1996a). The free- free 

component was fit assuming a spectral index (3ff = —2.15, the dust component was fit assuming 

(3d = +2.0. The coefficients given above combine the free-free and dust emission. 

c The coefficients give the most sensitive combination of the 31, 53 and 90 GHz data, in 

thermodynamic units, consistent with the constraint that emission with a spectral index (3/f = 

—2.15 (free- free emission) be nullified. The Haslam map is used to model synchrotron emission, as 

above. The DIRBE 140 fj,m map is used to model dust emission, as above. 
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Table 2. Scale-invariant Power Spectrum Normalization 



Map #l a 


Map #2 a Qrms-PS (/"K) x <2 


rms—PS 




9 

x 2 




oniric iiet^ueiicj ciuoo cuiieidLiuii 










31A 


31B 18.2 ± 4.1 68.3 


18.0 


± 


4.6 


71.0 


53A 


53B 18.3 ± 1.6 73.5 


18.6 


± 


1.7 


69.6 


90A 


90B 16.4 ± 2.2 72.3 


18.4 


± 


2.3 


71.0 




Single frequency auto-correlation 










31ws 


31ws 17.1 ± 3.7 67.9 


17.6 


± 


4.0 


79.6 


53ws 


53ws 18.7 ± 1.6 99.9 


19.4 


± 


1.6 


97.5 


90ws 


90ws 17.5 ± 2.0 63.5 


19.0 


± 


2.2 


61.5 




Multi- frequency cross-correlation 










53ws 


90ws 17.2 ± 1.5 60.8 


17.8 


± 


1.5 


64.1 


53ss 


90ss 17.0 ± 1.6 61.2 


17.9 


± 


1.6 


62.3 




Multi- frequency auto-correlation 










53+90 


53+90 18.5 ± 1.4 84.2 


19.6 


± 


1.5 


83.5 


31+53+90 


31+53+90 18.6 ± 1.4 80.0 


19.3 


± 


1.4 


78.2 




Multi-frequency auto-correlation with Galaxy 


model 






Correlation 


Correlation 17.5 ± 1.4 76.2 


18.5 ± 


1.4 


78.3 


Combination 


Combination 16.7 ± 2.0 89.2 


17.8 


± 


2.2 


92.4 



a The coefficients that comprise the map combinations in these columns are given 
in Table |I|, except for 53ss and 90ss which are straight sum maps: (A+B)/2. 

h (-min is the lowest order multipole remaining in the map after subtracting a best- 
fit multipole of order £ m in — 1- Qrms-PS is the most-likely quadrupole normalization, 
after calibrating the likelihood with Monte Carlo simulations, x 2 is tabulated, for 
reference, with respect to the mean of a scale-invariant model with the corresponding 
most-likely normalization. There are 71 bins in the 2-point function. 
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Fig. 1. — 2-point correlation functions obtained from the individual frequency maps, after monopole 
and dipole subtraction. The left-hand panels show the auto-correlation functions obtained from a 
weighted average of the A and B channel maps. The right-hand panels show the cross-correlation 
of the A and B channels, which are sensitive only to common structure in the maps. The error bars 
represent the uncertainty due to instrument noise, as described in the text. To guide the eye, the 
solid line is the auto-correlation of the weighted average of all six channels maps. All six 2-point 
functions are statistically consistent with each other. 
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Fig. 2. — 2-point correlation functions obtained from the multi-frequency maps after monopole and 
dipole subtraction, top left) The auto-correlation function of the weighted average map constructed 
from all six DMR channels, top right) The cross correlation function of the 53 GHz weighted sum 
map with the 90 GHz weighted sum map. bottom left) The auto-correlation function of the weighted 
average map with best-fit Galaxy template maps subtracted from the map (Kogut et al. 1996a). 
bottom right) The auto-correlation function of the linear combination map designed to cancel free- 
free emission. This map has a best-fit model of the synchrotron and dust emission also subtracted. 
In all panels the error bars represent the uncertainty due to instrument noise, as described in the 
text. To guide the eye, the solid line is the auto-correlation of the weighted average map (top left 
panel). 



